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Abstract. It is well known that the hydrodynamic limit of an interacting 
particle system satisfying a gradient condition (such as the zero-range process 
or the symmetric simple exclusion process) is given by a possibly non-linear 
parabolic equation and the equilibrium fluctuations from this limit are given 
by a generalized Ornstein-Uhlenbeck process. 

We prove that in the presence of a symmetric random environment, these 
scaling limits also hold for almost every choice of the random environment, with 
an homogenized diffusion coefficient that does not depend on the realization 
of the random environment. 



1. Introduction 

Consider a system of particles evolving on a multidimensional, periodic integer 
lattice of period 2N . Each particle performs a continuous-time random walk with 
rates p{x, y) that depend on both the position x and the destination site y. These 
rates are chosen as a fixed realization of a random field, in such a way that the 
resulting single-particle random walk is reversible with respect to the counting 
measure on the lattice. We call these rates the random environment. 

Particles interact between them only when they share a site, through an inter- 
action function g : No — » R+. The dynamics for this system is the following. At 
each time t, let rjt{x) denote the number of particles at the site x. For each pair 
of sites (x, y), after an exponential waiting time of rate g(n t (x))p(x, y) the particle 
at site x jumps to site y. This is done independently for each pair (x,y) and after 
each jump, the exponential waiting time for each pair (x, y) starts afresh. 

Such a system can be understood as a model for diffusion in heterogeneous 
media. The purpose of this article is to study the scaling limits of this system as 
N — > oo and mostly the influence of the randomness in this limit. As we will see, 
when the underlying random field is ergodic, stationary and satisfies an ellipticity 
condition, for any realization of the random environment the scaling limit depends 
on the randomness only through some constants which depend on the distribution 
of the random transition rates, but not on the particular realization of the random 
environment. 

In this article we study two related scaling limits for this process: the hydro- 
dynamic limit and the equilibrium fluctuations. The first one is a law of large 
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numbers for the empirical distribution of particles when the process starts from 
a configuration of particles with macroscopic density close to some initial profile 
while the second one is a central limit theorem for the empirical distribution of 
particles when the system starts from an equilibrium measure. 

The hydrodynamic limit has been obtained in the context of exclusion processes 
in [4] when the dimension d > 3 and extended in [14| to any dimension. In these 
references, it is not assumed the reversibility for the one-particle random walk with 
respect to the counting measure on the lattice, so in this sense their results are 
more general than ours. Their approach is based on the generalization of the non- 
gradient method of Varadhan [16] , [13] for the case of random transition rates. In 
[3], [8] the one-dimensional simple exclusion process is considered. 

In the reversible situation, we introduce the corrected empirical process. This 
process satisfies the gradient condition, which is a key property from which hydro- 
dynamics and equilibrium fluctuations can be easily obtained like in the non-random 
situation [6], [2]. Therefore, our approach is simpler, does not require any mixing 
condition and can be generalized to situations in which the non-gradient method 
does not apply, like kinetically constrained particle systems, the zero-range process 
with bounded interaction rate and particle systems in non-homogeneous lattices 

The introduction of a corrected empirical measure can be understood as a version 
of Tartar's compensated compactness lemma in the context of particle systems. In 
this reversible situation the averaging due to the dynamics and the inhomogeneities 
introduced by the random media factorize after introducing the corrected empirical 
process, in such a way that we can average them separatedly. For the dynamic av- 
eraging, we use the entropy method of [5] to derive the hydrodynamic limit, while 
for the equilibrium fluctuations we adopt Chang's proof of the Boltzmann-Gibbs 
principle [2] ; for the averaging of the random environment we use T- convergence. 
With this procedure, the scaling limits of the corrected empirical process are ob- 
tained. After this, we prove that in the limit as N — > oo, the corrected empirical 
process and the original empirical process are close enough to recover the scaling 
limit for the original empirical distribution of particles. 

In order to see how far can this picture be taken, we also prove the Boltzmann- 
Gibbs principle for functions that depend on both the particle configuration and 
the random environment. Notice that this more general version of the Boltzmann- 
Gibbs principle is not needed to obtain the equilibrium fluctuations for the empirical 
density of particles. 

The Boltzmann-Gibbs principle states that non-conserved quantities oscillate 
faster than conserved quantities, and therefore when averaged in time, only the 
projections over the density field are observed. In consequence, the Boltzmann- 
Gibbs principle is interesting by its own. In order to give further motivations for 
the study of the Boltzmann-Gibbs principle for random functions, we present two 
applications at the end of the article. 

The article is structured as follows. In Section [2] we describe the model and the 
main results. Section [3] is devoted to the proof of the hydrodynamic limit for this 
process and in the subsequent section we present the equilibrium fluctuations. The 
proof of the Boltzmann-Gibbs principle is referred to Section [S] For the reader's 
convenience, we include some well-known, but rather technical lemmas and defini- 
tions in the Appendix. 
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2. Notations and results 

2.1. The zero-range process. We define the zero-range process as a continuous- 
time Markov process rjt with state space ilff — {r/ : Tf^ — > No}, where Tfj is 
the cJ-dimensional discrete torus N~ lr L d /2N1 d . We consider as a subset of 
U d = [— 1,1] with periodic boundary conditions. This process has a generator 
whose action over local functions / : Qjy — » R is given by 

C N f(v)= E PN(x,y)g{ri(x))[f(v xy )-f(v)], 

where pn ■ ^% X — * K + is the jump rate of a random walk in Tff, g : No — » M+ 
is the interaction rate between the particles and j^" G flff is given by 

i7?(a;) — 1, z — x 
r)(y) + l,z = y . 
T){z),z ^x,y 

Notice that the dynamics of r\t conserves the number of particles. In particular, 
the process rj t is well defined for any initial configuration 770 G Ojv, since in that 
case the state space is finite. 

We will assume that the interaction rate g has linear growth: 

3c > : Cq 1 u < g(n) < c n Vn G N (2.1) 

We will also assume that the motion of a single particle is a nearest-neighbor 
random walk, so we take pjy(x, y) — if \x — y\ ^= l/N, where \x — y\ — J2i<d \ x i~ Vi\ 
is the sum norm in R d . This last hypothesis is not essential, but it simplifies the 
notation. We further assume that pn(x,v) — Pn(v,x) for all x,y G Tfj. This 
hypothesis will ensure the reversibility of the process r\t with respect to the measures 
v p defined below, and the reversibility of the randomm walk generated by pn(x, y), 
which is crucial in what follows. 

For each a > 0, let v a be the product measure in Qfy whose marginals are given 

by 

v a {v,v{x) = k} = ^r^r-^TT, 
Z{a) g{k)\ 

where g(k)l — g(l) ■ ■ ■ g(k) for k > 1, 17(0) = 1 and Z(a) is the normalizing constant 
for which i> a (^lff) = 1. By the linear growth of g (|2.1[) . v a is well defined for all 
a > 0. 

Define p = p(a) as the density of particles with respect to P a , namely: 

1 ka k _ aZ'(a) 



p{a) = E, a [ V (x)] = £ 



, , Q Z(a)g{k)\ Z(a) ' 

Again by the linear growth of g, a <— > p(a) is an homeomorphism from [0, 00) 
to [0,oo) and the inverse function a = a(p) is well defined for all p G [0, 00). We 
define v p = v a ( p ) and <p(p) — E„ p [g(r](0))]. Due to the symmetry of pjv(x,y), the 
measure v p is invariant and reversible for this process. 

2.2. The random environment. Now we discuss the choice of the jump rates 
pjy(x,y). Let (X^T^P) be a probability space and take a family {9 x \x G Z d } of 
^-"-measurable mappings 6 X : X — > X such that 

i) Pie- 1 A) = P(A) for all A G T, x e Z d . 
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ii) 6 Z Z , = 9 z+7j for all z, z' G 7L d . 

iii) If 9 Z A = A for all z G Z d , then P(A) = or 1. 

In this case we say that the family {0 x } X £Z d is invariant and ergodic under P. 
Let a = (oi, •••,d<i) : A" — > E d be a JF-measurable function such that there exists 
£o > with 

eo < a,i(u>) < €q 1 for all w G X and j = 1, d. (2-2) 

Fix lu G X. For each x G {-1 + 1/N,-1 + 2/N,..., l} d and i = 1, define 

p N {x,x + e l /N) =p N {x + ei/N,x) = N 2 ai {9 Nx u) (2.3) 

to which we call the random environment. 

For each G : Tfj — > K, define the operator LatG by 

L N G(x)= 53piv(aj,tf)[G(tf)-G(i)]. 
In the space of functions Zjv(Tjy) = {/ : — > R}, define the following norms: 

and 

ll/ll?,* = H/iiojv + ^ E ^[/(w) - /^)] 2 - 

|ar- v |=l/JV 

We denote by C 2 N the space of functions Zjv(Tjy) endowed with the norm || • | |o,jv 
and by (•, -)jv the inner product in C 2 N . Define Hi,n as the space of functions in 
InC^n) endowed with the norm || • 1 1 i,jv - 

Denote by C 2 (U d ) the space of square integrable functions in U d with respect to 
the Lebesgue measure and by || • ||q the corresponding norm in £ 2 (U d ). For each 
k > 1, denote by Hk(U d ) the Sobolev space in U d defined as the completion of 
Coo(U d ) under the norm 

n/ni= E mil 

|ce|<& 

where |a| denotes the order of the multi- index a and d a is the partial derivative of 
order a. 

The definition of convergence of a sequence Jn G Hi,n ( or £%) to / G TL\{U d ) 
(or C 2 {U d )) is given in Appendix IB. II 

From the homogenization theory the following holds: 

Proposition 2.1. Fix a typical realization o/pjv(v) and A > 0. There exists a 
positive defined matrix A that depends only on the distribution of a — (<2i,...,a<j) 
such that for any f^ and f such that /jv G H—i.n converges strongly to f G 
7i-i{U d ), un converges weakly in TLi.n to u, where u^ is defined as the solution 
of the equation 

\u N (x) - L N u N (x) = fn{x) 
and u is the solution of the equation 

Xu(x) - V • AS7u(x) = f(x). 
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A proof of this proposition can be found in [12] , Notice that the statement of 
this proposition makes sense for any choice of the jump rate pn{x, y). 

In order to prove the hydrodynamic limit we need this property on the jump 
rates pn{x,v) and for this reason we introduce the following definition. 

Definition 2.2. We say that a family of jump rates {pn ■ T^r xT^, — > K+jv admits 
homogenization, if there exist a constant eo > such that €qN 2 < pn{x, y) < e^N 2 
and a matrix A such that for any f G H.-i(U d ) smooth enough there exists a 
sequence /at converging strongly in H-i,n to f such that the solution Ujy G Hi,n 
of the equation 

Xu N {x) - L N u N {x) = f N {x) 
converges weakly in 7i\,N to the solution u G Hi(U d ) of 

Xu(x) - V • AVu(x) = f(x). 
In this case, we say that the matrix A is the T -limit of L^. 

For our purposes, / will be smooth enough if it is three times countinuously 
differentiable. 

Remark 2.3. By the theory of T -convergence, the matrix A satisfies the coercive- 
ness assumption eo|£| 2 < ^2ij £,i£,jAij < e^ 1 1^| 2 for all vectors £ G K d . In the 
previous definition, nothing excludes the possibility of the matrix A to be a function 
of the position x G U . See [llj for a one- dimensional example on which the T -limit 
of Ln is not constant in space. 

2.3. Hydrodynamic limit. Fix a function po : U d — > R+. A family of measures 
{pn}n>i in Q% is said to be associated to the profile po if for any function G G 
C{U d ) and any e > 0, 

lim fx N (rj G tl d N ; —[ v(x)G{x) - j p Q {x)G{x)dx > e) 

Here and in the sequel, denote by E PN the expectation with respect to /ijv and 
by K PN the expectation with respect to , the distribution of the process r]t 
starting from /ijy in V([0,T],Q%). We follow the evolution of the process rj t m a 
finite time interval [0, T] in order to avoid uninteresting complications due to the 
lack of compactness of [0, oo). 

Let p > be a fixed density. The entropy of p^ with respect to v p is defined by 



> e = 0. 



H N (p N \v p ) = 



I^og^ p ,ifp N «u p 
+oo, otherwise 

where for two measures p and v, p <C v means that the measure v is absolutely 
continuous with respect to p. 

We introduce a partial order < in £L d N as follows. For 77, 77' in we say that 
rj ^< 77' if rj(x) < r]'(x) for every x G T^. Once there is a partial order in the space 
state £1%, we can introduce a partial order in the space of measures in We say 
that pn is stochastically dominated by v p (also denoted by pjq ^ v p ) if there exists 
a measure p in flfj x flf^ such that: 

i) For all 77 G Q, N , p{v, &n) = MaK??)- 

ii) For all 77 G to N , p{^N,il) = v P (r))- 
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iii) The set {(i],rf))r) < 77'} has full measure under p. 
In this case we say that /I is a coupling of pn and Vp. 

Theorem 2.4. Let po : U d — > R be a bounded profile, and let {pn}n>i be a 
sequence of measures in ilfj associated to the profile po. Assume that the interaction 
rate g(-) is non- decreasing and has linear growth (see Section \2. Suppose that 
there exist constants Kq and p such that H^^vp) < K$N d and p^ ~<vp for every 
N large enough. Suppose also that the jump rates pj\f(x,y) admit homogenization 
with homogenized matrix A. 

Then, for every t < T , every continuous function G : U d — > R and every 5 > 0, 



lim P„ N 



jfd G(x)th(x) - I G(ii)p(l.„),lu\ - A 



= 0, 



where p(t, u) is the unique weak solution of the hydrodynamic equation 

id tP = \7 -(AV^p)) 

1 P(0, •)=*>(•)■ 



(2.4) 



In the sake of completeness we introduce the definition of weak solutions of 
equation (12.41) . 

Definition 2.5. Fix a bounded profile po : U d — > R. A bounded function p : [0,T] x 
U d — > R is a weak solution of equation \2.J$ if for every function G : [0, T] x U d — > R 
of class C^^T] x U d ), 

J J [p{s,u)d s G{s,u) + cf>(p(s,u))V ■ AVG{s,u)}duds 

+ [ p {u)G(0,u)du= [ p{T,u)G{T,u)du. (2.5) 
Ju d Ju d 

Let A4+ be the set of positive Radon measures in U d . The empirical measure 
71"^ is defined as the process in T>([0, T], M.+) given by 

?rf (du) = Vt(x)S x {du), 

where 5 X is the Dirac distribution at x. 

For G : U d -> R continuous, define 7rf (G) = J G(u)tt^ (du). The statement 
of Theorem 12.41 is equivalent to say that under P MJV the random variables 7Tj (G) 
converge in probability to J G(u)p(t, u)du for every G continuous and every t £ 
[0, T]. We will prove a stronger result for w^: 

Theorem 2.6. Under the hypothesis of Theorem \2.J\ ir^ converges in distribution 
in £>([0, T], M.+) to the trajectory p(t, u)du. 

Remark 2.7. Since p(t,u)du is a deterministic element o/£>([0, T], M.+), the con- 
vergence in distribution of irf implies its convergence in probability, from which 
Theorem \2.4\ follows. 
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2.4. Equilibrium fluctuations. Now we state a central limit theorem for the 
empirical measure, starting from an equilibrium measure v p . Fix p > and denote 
by S(U d ) the Schwartz space of infinitely differentiable functions in U d . 

Denote by y N the density fluctuation field, a linear functional acting on functions 
G e S{U d ) as 

{G) = E G(x){rit{x) - p). (2.6) 

xeT d N 

Notice that 

In this way we have defined a process in X>([0, T],S'(U d )), where S'(U d ) is the 
space of tempered distributions, which corresponds to the dual of the Schwartz 
space S(U d ). 

Theorem 2.8. Consider the fluctuation field y N defined above. Assume that the 
interaction rate g(-) has linear growth and that the jump rates admit homogenization 
with homogenized matrix A. 

Then, for every ti,...,tk 6 [0, T] and every Gi,...,Gk E S(U d ), the vector 
(y£[ (Ci)> ■■■■>yi[(Gk)) converges in distribution to (y tl (Gi), ...,y tk {G k )), where y t 
is the generalized Ornstein-Uhlenbeck process of characteristics <j)'{p)V ■ AS7 and 

3. Proof of Theorem 12.41 

By remark [2771 in order to prove Theorem 1 2. 41 it is enough to prove Theorem 12. 61 
The proof of Theorem 12.61 follows the standard lines of the proof of hydrodynamic 
limit by the entropy method for interacting particle systems. The route to proceed 
is the following: 

First we show that the distributions of irf in P([0,T], M.+) form a tight se- 
quence. Then we prove that the limit points of ir^ are concentrated on trajectories 
of measures absolutely continuous with respect to the Lebesgue measure in U d 
with a bounded density. Finally, we prove that these limit points are concentrated 
on weak solutions of the hydrodynamic equation (|2.4[) . By the uniqueness of these 
weak solutions on the space of bounded functions we conclude that tt^ has a unique 
limit point, concentrated on the trajectory with density p(t, u), where p(t, u) is the 
weak solution of equation (12. 4|) . Since the topology of convergence in distribution 
is metrizable, we conclude that the whole sequence ir^ converges to p(t,u)du. 

Unfortunately, this plan cannot be accomplished directly for ir^ , but for another 
auxiliary process, the corrected empirical measure, that we define below. 

Let A > be fixed. A function G : U d —> R is said to be regular if the function 
/jv € C 2 N defined by /v(a^) = XG(x) — V • AVG(x) converges strongly in H-i.n to 
AG — V • .4V. Notice that a sufficient condition for G to be regular, is G S C 3 (U d ), 
where C 3 (U d ) denotes the space of three times continuously differentiable functions 
on U d . 

Let G : U d -> R be regular. For each N > 1, define R x G(x) = AG(a;) - V • 
AVG(x) and G% : T d N -> R as the solution of 



XG x N {x) - L N G x N (x) = R x G(x). 



(3.1) 



8 



PATRICIA GONCALVES AND MILTON JARA 



By Lemma EU the following estimates hold: 

||G^||o,jv < A -1 1 |i? A C?| |o,at, (3.2) 

^ PN(x,y)[G*(y)-G x N (x)] 2 <\- x \\R x G\\l N (3.3) 

and 

||G^||oo,iv < A^H^GIU,^. (3.4) 
We define the corrected empirical measure 7r^ v,A by 



^' X ( G ) =4l Vt{x)G x N {x). 



N d 

Notice that ir^' X (G) is defined only for G regular, so 7T t ' A is not a well defined 
process in T>([0, T], M.+). Lemma \B . 1 1 shows that is a well defined process in 
the Sobolev space 7i~ k (U) for k > 3. However, this point will not be relevant for 
our proof of Theorem 12.61 

Since M.+ is separable and the vague topology in Wf+ is metrizable, in order to 
prove tightness of ir^ in £>([0, T},M+), it is enough to show tightness of ir^(G) in 
T>([0, T], K) for G in a dense subset of the set C(U d ) of continuous functions in t/ d . 
Therefore, it is enough to prove tightness of n^(G) for G regular. 

By Dynkin's formula, 

M?(G) = ^ A (G) - ttq' x (G) - J* i Yl 9(vs(x))L N G x N (x)d s (3.5) 

is a martingale of quadratic variation given by 

(M?(G)) = J* j^s E 9(Vs(x)) P N(x,y)[G x N (y)-G x N (x)] 2 ds. 

x,yeTf, 

We claim that Mf (G) goes to as TV -> oo in £ 2 (P Miv ). In fact, 
E MJV [.Mf(G) 2 ] =E^[(^f(G))] 

= / j^d L E m [g( Vs (x))}p N (x,y)[G x N (y)-G x N (x)} 2 ds 

^^A-H^GH^^O. 

In order to obtain this last bound, we have used the estimate (|3.3p . the fact that 
jUy is stochastically dominated by Vp and Proposition lA.il 

To prove tightness for the martingale M.^(G), we use the following criterion, 
due to Aldous: 

Proposition 3.1. A sequence of probability measures {Pn}n in f([0, T],R) is 
tight if 

(i) For all < t < T and for all e > there exists a finite constant A such 
that supjv Pj\r(|#t| > ^4) < e, 
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(ii) For alle>0, 

lim limsup sup Pn{\x t +p ~ x T > 5) = , 

[3<S 

where T is the set of stopping times with respect to the canonical filtration, 
bounded by T . 

A proof of this lemma can be found in [S]. In our case, condition i) follows from 
the fact that M^(G) converges to in £ 2 (P MJV ) and Tchebyshev's inequality. On 
the other hand, by Doob's optimal sampling theorem, we have that 



?*s[\M?+i3(G)-M?(G) \ >e] < 



<^EW 



(M? +fl (G))-(M?(G)) 

T+/3 1 



E 



g(r, s (x))p N (x, y) [G x N (y) ~ G x N (x)] 2 da 



(3C(G,c ,\,e ) 
- ]yd-2 i^n 



E Vix) 



In this last bound we have used the conservation of the number of particles, the 
estimate (13.31) . and the uniform bound for p^{x,y). Since the expected initial 
density of particles is bounded by p, condition ii) follows. 
Notice that the integral term in (|3.5[) can be written as 



1 

N 1 



£ g{vs(x))[XG x N (x)-R x G(x)]ds. 



(3.6) 



We see that 



sup 

|s-t|<<5 



i g(vAx))[XG x N (x) - R x G(x)]dt' 



<5C(G,g)E^ N [ sup 
te[o,T] 



< 



\m\o,N > 



that goes to as 6 — » 0, uniformly in N by Lemma fA.21 Therefore, by Arzela-Ascoli 
criterion, the integral terms in (|3.5p form a tight sequence in D([0, T],R) and their 
limit points are concentrated on continuous trajectories. By equation (|3 . 5|) the se- 



quence tt^' X (G) is tight in Z?([0, T], R). On the other hand, since M.f (G) goes to 
in £ 2 (P MJV ), any limit point of A4^(G) has null finite-dimensional distributions. 
Therefore, Ai^(G) converges to in distribution as a process in 2?([0,T],R). Con- 

(G) are concentrated on continuous trajectories. 



sequently, the limit points of n t 
Notice now that 



i m [ sup 

te[o,T] 



n?' x (G)-n?(G)n <\\G 



G||o,v e m«[ SU P II^IIo.jv] 



~te[o,T] 

By Proposition ^. 11 ||G^ — G| |o,jv converges to as N — » oo, and by Lemma fA.21 
E^ N [sup f ||»?t||ojv] i s bounded in N. Therefore, sup t |7rf' A (G) — ir^(G)\ — » in 
£ 2 (P MJV ). A simple e/3 argument allows us to obtain from this result that 7Tj (G) 
is also tight in 2?([0,T],R) and that wj* (G) and 7r^' X have the same limit points. 



10 PATRICIA GONCALVES AND MILTON JARA 

Since the set of regular functions is dense in C(U d ), this ends the proof of tightness 
for Trf in V([0,T],M+). 

Let 7Tt be a limit point of irf , and let Q be its distribution in T>([0, T],Ai + ). For 
any positive function G 6 C(U d ), 

Q(w t (G) > M) < liminf Q(ir?(G) > M) 

N^OQ 



= liminf fj, N {N- d ^ rj t (x)G(x) > M) 
xeT% 

< \imM Vp(N- d v(x)G(x) > M) 



N— >oo 

< 1 (J G{u)du > M/p 

Here we have used once more, the fact of p^ being stochastically dominated by 
an invariant measure Vp and Proposition I A. ll 

Therefore, if < G < 1 then Q{ir t {G) > 2 d p) = 0. By the dominated conver- 
gence theorem, for every closed B C U d it holds that Q(n t (B) > pA(B)) = 0, where 
A denotes the Lebesgue measure in U d . In particular, the process ir t is concentrated 
on measures absolutely continuous with respect to A. 

Let 7r(t, u) be the density of 7r t with respect to A. The same estimates prove 
that Ti(t,u) is bounded by p in [0,T] x U d . 

Notice that R X G is a smooth function, but XG x f (x) it is not smooth. However, 



N 

± £ 9(Vs(x))[G x N (x)-G(x)]ds 

xtTt 



< 



ds 



: - 2 t f E MJV [(-^ r) s (x)\G x N (x)-G(x) 
<Co 2 t 2 [QzY, ri{x)\G x N {x)-G{x)\) 2 dv- p {ri) 

<c^J v (ord^\\G x N -G\\l N ^^0. 

In the previous we used Schwarz inequality together with the translation invari- 
ance of Vp. As a consequence, 



^f(G)=7rf' A (G)-7r ^ A (G)- 



/ 4d E 9{Vs(x))V ■ AVG(x)ds (3.7) 



plus a rest vanishing in £ 2 (P MN ) as N — + oo. The next result will allow us to write 
the integral term (|3.6p as a function of tt^ plus a vanishing term as N — > oo. 

Proposition 3.2. (Replacement Lemma) 



For every 8 > 0, 



lim sup lim sup I 



W E V eN{Vs,x)ds > S 

•JO 



xeTt 



= 0, 
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where 

V t (rj,x)= - ^2 g(r](x + y)) - <t>(rf(x)) 



(21-, ... 

\y\<i 



v 7 i»i<i 



The proof of this proposition is the same as the one presented in Chapter 5 of 
[9j, so we omit it. Using this proposition, we see that for any continuous function 

G : U d -> R, 

in P MJV -probability as N — > oo and then e — ► 0. On the other hand, since rf s (x) — 
tt^(1(|u — x\ < e)), we conclude that 

J -j^d 9(Vs(x))G(x)ds -> J ds J <t>{-K{s,u))G{u)du 

in P MN -probability. Since Aif(G) converges to 0, taking N — -> oo in equation (|3.7[) 
we obtain that 

= J n(t,u)G(u)du- J po(u)G(u)du - J J ^(tt(s, u))V • -4VG(u)dwds (3.8) 

for every G regular. Approximating a twice-differentiable function G by regular 
functions G n in the uniform topology, we extend this identity to functions G G 
C 2 (U d ). 

Let G : [0,T] x U d -> R be of class C 1 - 2 . Take the partition {i* = Ti/n;z = 
0, n} of the interval [0, T] and define G„ : [0, T] x U d -> R by 

G„(t,u) = G(*i_i,u) + G(ti,u), 

for i G [U-i, U]- In general, for a piecewise-differentiable path G : [0, T] — > Z^-, 

Trf (G t ) - Trf (Go) - f {*?(d.G a ) + ^ £ g{Vs(x))L N G s (x)}ds 

is a martingale of quadratic variation 
pt i 

/ X! PN(x,y)g(v s (x)) [G s (y) - G s (x)] ds. 

Repeating the arguments in the proof of equation (|3.8p for G n , we conclude that 
= f Tr(t,u)G n (t,u)du — j po(u)Gn(0,u)du 

{ir(s, u)d s G n (s, u) + 4>(n(s, it))V • AVG n (s, u)duds. 
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Taking the limit as n goes to oo, we obtain that 
= / n(t,u)G(t,u)du- / p o {u)G(0,u)du- 



o 



{n(s, u)d s G(s, u) + 0(7r(s, u)) V • AVG(s, u)duds 



for every G : [0,T] x C/ d — > R of class C 1,2 . This is the weak form of the hydrody- 
namic equation (|2.4p , see (12. 5p . Since equation (|2.4p has at most one weak solution, 
we conclude that 7r(i, w) = p(t, u) Q — a.s., which ends the proof of Theorem 12.61 



4. Proof of Theorem 12.81 

Denote by Q N the distribution in V([0, T],S'(U d )) induced by the process y t N 
and Pp. The standard proof of equilibrium fluctuations cannot be accomplished for 
the density field y N . In order to overcome this problem we introduce as before, 
the corrected density fluctuation field defined on functions G £ S(U d ) by 

where is the solution of equation (|3.ip . 

For t > 0, let T t be the cr-algebra on D([0, T], S'(U d )) generated by y s (H) for 
s < t and H in S{U d ) and set T — cr(Ut>o Denote by the distribution on 
V([0,T],S'(U d )) induced by the corrected density fluctuation field y N ' X and v p . 

We make use of the following result, which permits to identify the limiting pro- 
cess: 

Proposition 4.1. 

There exists a unique process yt in C([0, T},S'(U d j) such that: 
i) For every function G £ S(U d ), 

M t (G) = y t (G) - y (G) - f 34(<^'(p)V • AWG)ds 



(M t (G)y ~ <j>( P )t / VG(«) • AVG(u)du 
Ju d 

are J- t -martingales. 
ii) y$ is a Gaussian field of mean zero and covariance given by 

E[y (G)y (H)] = X (p) [ G{u)H{u)du, (4.1) 
Ju d 

where x{p) = Var(n(Q), v p ) and G, H £ S(U d ). The process yt is called the 
generalized Ornstein-Ulenbeck process of mean zero and characteristics 0'(p)V-,4V ; 



Theorem l2.8l is a consequence of the following result about the corrected fluctu- 
ation field. 

Theorem 4.2. Let Q be the probability measure on C([0, T], S'(U d )) corresponding 
to the stationary generalized Ornstein- Uhlenbeck process of mean zero and charac- 
teristics <f)'(p)V ■ AV, y/ (f>(p)AV. Then the sequence {Q^}n>i converges weakly 
to the probability measure Q. 
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Before we enter into the proof of this theorem, we prove Theorem 12.81 from it. 
In fact, it is enough to show that 

lim E„, [(y t N (G) - y?' X (G)) 2 ] = (4.2) 

iv — > oo 

for any t G [0,T], G G S(U d ). But this is immediate from the fact that 
converges to G in C 2 N and the independence of r)(x), r)(y) for x ^ y under the 
invariant measure v p . 

In order to prove Thcorcm l4.21 we need to verify that the sequence of probability 
measures {<2^}jv>i is tight and to characterize the limit field. Then we show that 
the limit field is equal in distribution to y t using its characterization in terms of 
the martingale problem (Proposition 14. 1] ) . 

Fix a smooth function G G S{U d ). By Dynkin's formula, 

M t N ' x (G) = y t N '\G) - y?>\G) -j'jl^Yl 9{ris(x))L N GUx)ds (4.3) 

is a martingale with respect to the natural filtration Tt — a(j] s ,s < t) whose 
quadratic variation is given by 

(M t N ' X (G)) = f Y, 9{r, s {x))p N {x,y)[G%{y)-G x N {x)] 2 ds. 

At first, we establish the limit of the quadratic variation. Notice that in the 
previous formula we can replace g\r) s {x)j by <j>(p), since 



± Y {9(Vs(x))-cl ) (p)}p N (x,y)[G x N (y)-G x N (x)} 2 d S 



< 



- ATd~2 II- rAg 'IIo,vII^ a IIoo,v- 



< J^VarQ?,^) Y PN(x,y)[G x N (y)~G x N (x)] 2 x 

x sup Y PN{x,y)[G x N {y) - G x N (x)] 2 

xeT " y&K 
Ct 2 

]\[d-2 

For dimension d > 3, this last expression goes to as N — > oo. In order to cover 
the case d = 2, we can use Theorem (1.31) of [15] . expression (1.32) with t — s and 
a = 1/N and take the Laplace transform of equation (1.32), to obtain a sharper 
estimate for G^(x) — G^(y). In this case, we obtain that the last line is bounded 
by j\f-( d - 2 + 2a ) ; for some a > 0. As a consequence, for any d > 2, the quadratic 
variation can be written as 

/ ^ E <Kp> E PN(x,y)[G x N (y)-G x N (x)] 2 ds, 

x£T d N y£T d N 

plus a vanishing term in the £ 2 (P„ )-norm. Using the convergence of G^ in C 2 N 
and the resolvent estimates in the proof of Lemma \B.1\ this last integral converges 
to 



t<j)(p) [ VG(u) • AS7G(u)du, 
Ju 
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as N goes to oo. 

Now we study the limit of the martingale M t Ar ' A (G), see expression (|4.3p . Since 
Y^xeT d I j nG^ 1 {x) — 0, we can rewrite the integral part of the martingale as 

fil MVs(x)) - <p{p)}L N G x N {x)ds. 

On the other hand, since G^ is the solution of equation (|3.ip . the last integral 
can be written as 

<■* 1 „ . r . y, 

>ds. 



N d/2 

xert 



{g(Vs(x)) - <f>( P )}{\G x N (x) - XG(x) + V • AVGix)^ 



Our aim now consists in showing that it is possible to write the integral part 
of the martingale as the integral of a function of the density fluctuation field plus 
a term that goes to zero in £ 2 (P !/ ). The first result needed to proceed in that 
direction is the following: 



1 



o N d / 2 



{g(vs(x)) - Hp)}[g x n (x) - G(x)]ds 



< 



xETt 



<CVar( g) u p )\\G x N -G\\l N 



N- 



The second one is known as the Boltzmann-Gibbs principle. Here we have the 

X. For each 



N 



need to introduce some definitions. Take a function / : \ x flf^ 
uj G x an< i eacn x G T^, define 

f(x, rj) = f(x, t), (j) =: f{9 xN u, t x t}), 

where t x t) is the shift of rj to x: T x rj(y) = rj(x + y). Notice that we do not include 
explicitly the dependence of f(x,r]) in u), since in our setting u> is fixed. 

Definition 4.3. We say that f is local if there exists R > such that f(u),ij) 
depends only on the values of rj(y) for \y\ < R. In this case, we can consider f as 
defined in all the spaces \ x f or N > R. 

Definition 4.4. We say that f is Lipschitz if there exists c = c(u>) > such that 
for all x, \f(u>,r)) — f(u>,ri')\ < c\rj(x) — Tj'(x)\ for any r\, rj' such that r/(y) = r/'(y) 
for any y ^ x. If the constant c can be chosen independently of u, we say that f is 
uniformly Lipschitz. 

Theorem 4.5. (Boltzmann-Gibbs principle) 

For every G G S(U d ), every t > and every local, uniformly Lipschitz function 

f: X xn d N ^m, 



lim 



N d/2 



G(x)V f (x,r) s )ds 



= 



(4.4) 



where 



Vf(x,rj) = f(x,r)) - E Vp [f(x,r))] -d p E j f{x,rj)dv p (ri) (p(x)-p). 
Here E denotes the expectation with respect to P, the random environment. 
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In order to simplify the exposition, the proof of this last result is postponed to 



the next section. As we need to write the integral part of the martingale M t : (G) 
in terms of the density fluctuation field, by using the first result stated above we 
are able to write the integral part of the martingale as 



1 



N d/2 



{g{r} s [x))-<t>{p)}V-AVG{x)ds 



plus a term that converges to in the C 2 {T V )-norm. The replacement of the 
function <7(?7 S ) — <fi(p) by </>' (p)[rj s (x) — a] in the last integral, is possible thanks to 
the Boltzmann-Gibbs principle. Doing so, the integral part of the martingale can 
be written as 

\G)-y^\G)~ JL £ <l>>(p)V-AVG(x)(r) s (x)-p)d S 



N.X, 



plus a term that vanishes in £ 2 (P J/ ) as N — > oo. Notice that the integrand in the 



previous expression is a function of the density fluctuation field yf , see (|2.6I) . By 
(|4.2p . we can replace inside the integral of last expression the density fluctuation 



field by the corrected density fluctuation field y t ' . 

Suppose that the sequence {Q^}n>i is tight and let Q x be a limit point of it. 
Denote by 3^t the process in T>([0,T],S'(U d )) induced by the canonical projections 
under Q x . Taking the limit as N — » oo under an appropriate subsequence in 
expression (]4.3|) . we obtain that 



M t \G) = y t (G) - y (G) 

is a martingale of quadratic variation 



y s {(t>\p)V ■ AVG)ds 



t<p(p) / VG(ti) • AVG{u)du. 



On the other hand, it is not hard to show that 3^o is a Gaussian field with 
covariance given by (|4.ip . Therefore, Q A is equal to the probability distribution Q 
of a generalized Ornstein-Uhlenbeck process in C([0, T], S'(U d )) (and it does not 
depend on A). 



As a consequence, the sequence {Q^}n>i has at most one limit 



point and Theorem 14.21 shall follow if we prove tightness for {Q^}n>i 

Lastly, it remains to treat the problem of tightness of the sequence {Qjy} 



JV>1- 



For that we use a criterion due to Mitoma [TD] (see also [S]), which allows to conclude 
that the sequence is tight and that any weak limit is supported in C([0, T], S'(U d )), 
since the following estimates hold: 

a) For every T > and G £ S{U d ), 



sup sup J 

N te[0,T] 



< oo. 



sup sup J 

N te[0,T] 



N d/2 



4>'{p)V ■ AVG(x)(r) s (x) - p) 



< OO. 



sup sup J 
n te[o,T] 



1 



Nd/2 E 



ftp) PN(x,y)[G x N (y)-G x N (x) 



< oo. 
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b) For every G G S{U d ) there exists S(t, G, N) such that linijvwoo S(t, G, N) 
and 



liml 

N 



Vp { sup y^ x (G) - y N s : x (G) > S(t, G, N)) = 



The first expectation in a) is bounded by | |G^| \ooX(p) i which in turn is bounded 
by C||i? A G|| 00 . The second expectation in a) is bounded by C||V ■ AVG\\% and the 
last one bounded by C\\R X G\\%. 

To prove b) we only have to remark that by definition of the process it holds 



that su Po < s < t |3^ A (G) - y^(G) 

By the results proved {Q^f} n >i is tight and we have identified above a unique 
limit point Q that corresponds to the Ornstein-Uhlenbeck process; consequently 
the whole sequence converges to Q. 



< 



N d / 2 



5. BOLTZMANN-GlBBS PRINCIPLE 

This section is devoted to the proof of Theorem (|4.5p . Let / : \ x ^jv — * +°° 
be a local, uniformly Lipschitz function and take f(x,rf) = f(0NxU, T x rj). 

Fix a function G G S(U d ) and an integer K that shall increase to oo after N. 
For each N, we subdivide in non overlapping cubes of linear size K . Denote 
them by 1 < j < M d }, where M — [t^]. Let Iq be the set of points that are 
not included in any Ij which implies that \Ij\ < dKN^ 1 . If we restrict the sum 
in the expression that appears inside the integral in (|4.4j) to the set Iq, then its 
L 2 (P„ ) norm clearly vanishes as N — > +oo. 

Let A s/ be the smallest cube centered at the origin that contains the support of 
/ and define Sf as the radius of A Sf . Denote by I® the interior of the interval Ij , 
namely the sites x in Ij that are at a distance at least s / from the boundary: 

I° = {x€lj,d(x,T d N \I j )>s f }, 

Denote also by I c the set of points that are not included in any By construc- 
tion it is easy to see that \I C \ < dN d (^- + Using the notation just settled, we 
have that 



N d ' 



x — 2 J2 mv f (x,vt) = ^E 



x,Vt) 



x£l c 



M a M a 

—2 E E [ h w - H ^)] v f^ + wi* E E v t& 



N d / 



where yj is a point in Ij. We assume that the points yj have the same relative 
position on each of the cubes. The first step is to prove that 



lim lim Ej, 

K — >oc A'' — >oo 



1 



N d/2 



J2 H(x)V f (x,r, t )ds 



xei c 



0. 



Applying Schwarz inequality, since v p is an invariant product measure and since 
Vf has mean zero with respect to the measure v p , the last expectation is bounded 
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above by 



]T H(x)H(y)E Up [V f (x,r ] )V f {y,r ] )]. 



x,y£l c 
\x-y\<2s f 

Since V f belongs to C 2 {v p ) and \I C \ < dN d { K 
ishes by taking first N — > +00 and then K — > +00. 

Applying the same arguments, it is not hard to show that 



d/£(Jl -J- the last expression van- 



lim E„ 

N^oc 



1 



N d / 2 



E E [^)-^(%)]V>(a:,r 7t )d S 



3=1 

In order to finish the proof it remains to show that 



lim lim E,, 

K — »oo N — >oc 



1 



A I 



N d/2 



E^(%) E v fMds 



3 = 1 



= 0. 



Let L^v be the generator of the zero-range process without the random environ- 
ment (that is, taking a(uj) = 1 in (|2.3[1 ). and without the diffusive scaling iV 2 . For 
each j — 0,..,M d denote by Q the configuration {rj(x),x € Ij} and by hi, the 
restriction of the generator L^r to the interval Ij, namely: 



1^)= E 9(v(x))[h(r 



x,yelj 
\x—y\ = l/N 



We point out here that we are introducing a slightly different generator than the 
one that generates the dynamics, namely Cn- The reason for doing this stands on 
the fact that the dynamics generated by this operator is translation invariant. The 
generator that we choose to introduce here is not random, but due to the ellipticity 
assumption on the environment, it is mutually bounded with the one that we have 
started with. 

Now we introduce some notation. Fix a local function h : \ x ^v — ¥ ^> mea- 
surable with respect to a(i](x),x G I\), such that E[J h(u),Tii) 2 dvp\ < 00 and let hj 
be the translation of h by yj — yo: hj(x,rj) = h(9( yj - yo )NUj,Ty j ^ ya r/). Denote by 
C 2 (i>p x P) the set of such functions. Consider 

A I ' 



3 = 1 

By proposition A 1.6.1 of [5] and the ellipticity assumption, it is not hard to 
show that 



E„ 



1 



AI" 



o N d / 2 



<20 eo -4|||^||| 2 _ 1 , 



3=1 



where the norm 



is given by the variational formula 



\V N III 2 
I v H.h\\\-l 



sup 



{2 f v£ h (r,)F(r))dv p - N 2 (F, -h N F) p }, (5.1) 



where (•, ■) p denotes the inner product in L 2 (v p ). 
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By the Cauchy-Schwarz inequality, 

UMQ F (v)dv P < ■±p{-'Li j h j ,h j ) p + ^(F, -U 3 F) P 

for each j, where jj is a positive constant. Therefore, 

f 2 Md 1 

2 / Vi th (v)F( V )du p < —^H(y j ){ — (--L Ij h j ,h j ) p + ^(F i -h Ij F) p }. 

3=1 3 

Taking for each j, jj — N 2+ i\H{yj)\^ 1 we have that 
i M " , 



jyd/2 

3=1 

and the expectation becomes bounded by 



1, M 



E ^(y^i . , , " ii" n-x 1 t ; , . 

— ( — Ljj ftj > fyj )p < e2N2+d 2^m*( j j ' j) 



20eW x t 



3=1 ,J 



20iAf d ||i?| 



3=1 



By the ergodic theorem, the sum in the previous expression converges as N — > oo 
to a finite value and therefore this last expression vanishes as N — > oo. To conclude 
the proof of the theorem we need to show that 

lim inf 

K—toc heC 2 (y p xP) 



lim E„ 



ft 1 M 2 

3=1 KG/? 



3=1 xe^ 

By Schwarz inequality the expectation in the previous expression is bounded by 

e M " 



3 = 1 



167° 



because the measure ^ p is invariant under the dynamics and also translation invari- 
ant and the supports of Vf(x,rj) — Lj 4 /i*(&) and Vf(y,i]) — hj.hj(Q) are disjoint 
for x £ J° and y S 1°, with i 7^ j. 

By the ergodic theorem, as N — > 00 this expression converges to 



K d 



\H\\tE 



(5.2) 



So it remains to show that 



inf S 

de£ 2 (i/ p xP) 



( E ^/^'^ - ^iM^,v) 



dv, 



0. 



Denote by T^QL/j) the range of the generator L/j in C?(y p x P) and by /^(Lj-J 
the space orthogonal to TZih^). The infimum of (|5.2[) over all ft- G C?(v p x P) is 
equal to the projection of ^2 x£l o Vf(x,rf) into TZ^h^) . 
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It is not hard to show that TS-QLtJ is the space of functions that depends 
on r\ only through the total number of particles on the box I\. So, the previous 
expression is equal to 



K^oo K d 



n 



dv c , 



(5.3) 



where 7] h = K d Y, x eh ^i?)- 

Let us call this last expression 2"o- Define ip(x,p) — E Up [f(Q x Lo)]. Notice that 
Vf{x,n) — f(x,ri) — ip(x,p) — E[ip'(x, p)] (n(x) — p), since in the last term the 
derivative commutes with the expectation with respect to the random environment. 
In order to estimate the expression (|5.3|) using the elementary inequality (x + y) 2 < 
2x 2 + 2y 2 , we split it into three pieces: To < 4(Xi + T 2 + ^3), where 



X\ = —;E 



Y,E Vp [f{x,ri)W\ 

x£l° 



ipix.'q 1 



dv r 



E 



la 



K d 



E 



x£l° 



We will make use of the following lemma, known as the equivalence of ensembles. 

Lemma 5.1. Let h : flfj — » R a local, uniformly Lipschitz function. Then, for each 
[3 > there exists a constant C that depends on h only through its support and its 
Lipschitz constant, such that 

C 



E v \h{n)\if]-E„ N [h{rj)\ 



< 



N d 



whenever p, 77 < (3, where 



V 



N 



\x\<N 



In order to estimate I\ and T 2 , we introduce the indicator functions l(r/ /l < (3). 
By a large deviations estimate, v p (j] Ix > (3) < exp(—C(f3)K d ). Since / is Lipschitz, 
it has bounded exponential moments of any order and a simple Schwarz estimate 
shows that we can introduce the indicator function 1(?7 71 < (3) into the integrals in 
X\ and Xi . By Lemma 15.11 



1 



E 



x€l° 



< 



c_ 
K d- 



which vanishes as K — > 00. 

Using a Taylor expansion for ip(x,p), we see that 

1 



K d 



E 



and also goes to as K — > 00. 



< 



C_ 

K 1 
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Finally, we see that 

2 3 = E Vp [(r,(0) - p) 2 } -E[{±- d P) ~ EW{x, p)] 

X 

and it goes to as K — > oo by the £ 2 -ergodic theorem. 

5.1. Some applications of the Boltzmann-Gibbs principle. In the proof of 
Theorem 12.81 we need to use the Boltzmann-Gibbs Principle 14.51 for the function 
17(77(0)), that does not depend on the random environment. In particular, the results 
of the previous section are not needed in the proof of Theorem I2.8[ since the proof 
for the non-random case applies directly for functions that do not depend on the 
random environment. We point out here two applications for the Boltzmann-Gibbs 
principle as stated in Theorem 14.51 

First application: Consider, for simplicity, some local, bounded and uniformly 
Lipschitz function f(w,rj) that does not depend on the value of t](0). For each 
rj £ Ojv, define 

rj{x) + 1, z = x 
rj(z),z ^ x, 



r](x) — 1, z = x 
J](z), z ± x. 

Notice that 0~?7 is well defined only if rj(x) > 1. We can define a reaction-diffusion 
model adding to the zero-range dynamics a Glauber dynamics as follows: 

C r N d F( V )=:N 2 C N + f(x,v)[F(e+r,) - F( V )} 



Z^t^F'"- 1 '' F ""-' 



where we define f(x, Q^i])/g(ri(x)) — if rj(x) = 0. We have chosen the annihilation 
rate in such a way that the measure v p is invariant for this process. Therefore, we 
can obtain the equilibrium fluctuations for this model as in Section 2] 

Second application: This one has to do with the convergence of additive 
functionals of Markov processes. For each / satisfying the conditions of Theorem 
14.51 define the density fluctuation field for / acting on functions G € S(U d ) as 

Z t J (°) = ^ E G{x)){f(x,r] s ) - E Vp [f(x,r})]}. 

Note that for f(x, rj) = r](x) — p, the density fluctuation field for / is the density 
fluctuation field introduced above and denoted by yf(G). 
For fixed / as above, define the additive functional 

lf{t) = j^Z»{G)ds. 

Then, by Theorems dE] and EM 

lim Xf {t) = d p E I f(u>,r/)dis p / y s (G)ds in distribution. 

N^oc I J 1 J 
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Appendix A. Some estimates for r\ t 

A.l. Entropy production. Denote by /ijv(i) = [in the distribution of r\t in 
flff under P MJV and define / t w = . The density satisfies the Kolmogorov 

equation 

d 

~dt" 



-J t N (r ) )=£ N f t N (r 1 ). 



For each density / : flff — > R + , define the Dirichlet form T)^(f) by 
ZM/)= E J 9Hx))WJW)-VJWdvp, 

x,y£T N 
\x-y\ = l/N 

and the entropy Hjsr(f) = J f\ogfdvp. By the ellipticity assumption in ppj(x,y), 
the entropy production is bounded by the Dirichlet form of / t [9] : 

±H N (f t N )<-2e N 2 V N (f t N ). 

Assume that H N (fi N \^p) < K Q N d , or in other words that H N {f^) < K N d . 
Since the Dirichlet form and the entropy are convex functions of /, integrating 
the previous inequality we obtain the bounds 



where 



A. 2. Attractiveness of r/ t . Take two probability measures fx, v in f2jv such that 
[i < v. When the jump rate g(-) is non-decreasing, it is possible to construct a 
process (rjt,rj' t ) in fijv X fijV) starting from a coupling p, of /i and z/, such that for 
every t € [0, T] 

i) The distribution of r\ t in £>([0, T], S7jv) is equal to P„. 

ii) The distribution of r)' t in £>([0, T], fijv) is equal to P„. 

iii) The distribution of (r)t,rj' t ) in 2?([0, T], f2jv x f2jv) is concentrated on the set 
{(77,77') en N x fl N ;n < 77'}. 

In this case the process r] t is said to be attractive. We say that a function 
h : fi^v ~> K is non-decreasing if for 77 ^ 77' then /i(r/) < h{rf). The following 
proposition is an immediate consequence of the existence of the process (j] t ,r]' t ). 

Proposition A.l. Let p, v be two probability measures in fi^ such that fi ^ v. 
Let h : Q% -> M be a non- decreasing function. Then, 



for all t € [0,T]. 



E„ [/»(%)] < E„[Ji(»fc)] 
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A. 3. An C 2 estimate for rjt. Consider the process rjt starting from the equilibrium 
measure v p . Define the /^-norm of rjt by 



\Vt\\l, N 



Jd E r ^ a 



N d 



By Dynkin's formula, 



M 



N 



\Vt\\l,N ~ \\Vo\\o,N ~ / £N\\Vs\\l,N ds 



is a martingale of quadratic variation 



(A.l) 



(Mf)= / IM|^v(-£aOIMI^. 
Jo 

Explicit computations show that E l/p [(A4f r )] < C/N d ~ 2 . Therefore, by Doob's 
inequality, 



E„ 



sup \M t 
te[o,T] 



AT|2 



< C/N 



d-2 



For the integral term in !A.l( we have the following estimate: 



sup / C N \\r] s \\l N dsY 
' te[o,r] Jo 



< CtE v 



\H\IA~^)M\In 



Therefore, for dimension d > 2, we conclude that E„ p [sup t ||??t|lo at] is uniformly 
bounded in N. Since ||?y||g N is an increasing function, we have proved the following 
result: 

Lemma A. 2. Fix p > 0. Let {/^n}n>i be a sequence of measures such that 
fiN ^ v p for all N . Then, 

supE pjv [ sup ||%||o,jv] < +°°- 

N&~i t£[0,T] 



Appendix B. Functional analysis in the spaces C 2 n , Hi,n 

B.f. Convergence in C 2 N , Hi,jv. 

Fix / G Hi n- We define the linear interpolation Thf of / as follows. To 
fix ideas, take d = 3. We divide each of the cubes of size 1/N in Tfj into six 
tetrahedrons with vertices in Tjy. The way we do this is not important, but we do 
it in the same way for every cube in T^,. 

For a point u in one of such tetrahedrons, we define T^f(u) as the linear inter- 
polation of the values of / on the vertices of the tetrahedron. In this way we have 
defined a function T^f in Hi(U d ). 

We say that /a/ G Hi y N converges strongly (resp. weakly) in Hi. at to / G H.i(U d ) 

if 

lim T^fjq = f strongly (resp. weakly) in Hi(U d ). 

N— >oo 

In an analogous way, for each u G U d we define T$f(u) — f(x) if \u — x\ < 1/2N. 
We say that /at converges strongly (resp. weakly) in C 2 N to / G C 2 (U d ) if T^Jn 
converges strongly (resp. weakly) to / in C 2 {U d ). 
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A sequence Jn € H-x,n converges to / 6 H-i(U d ) strongly (resp. weakly) if 
for any sequence gM £ Ti-i,N and g 6 Tii(U d ) such that gM — > weakly in (resp. 
strongly) Wi,jv we have 

lim (/jvj 9n)n = (f,g)- 

B.2. Resolvent estimates. 

Let / be a regular function and let un be the solution of the resolvent equation 

Xu N (x) - L N u N (x) = f(x). (B.l) 

Lemma B.l. There exists a constant c = c(A) such that 

max{||uAr||o : Ar, ||uiv||i,jv} < c||/||i. 

Proof. By Lax-Milgram's lemma, this equation has a unique solution in Hi,j\r. Tak- 
ing the inner product of equation (|B.1[) with respect to un, we see that 

A||wjv||o,at + ( u n, ~L N u N ) N < (/, u N ) N . 

By the Cauchy-Schwarz inequality, |(/, itjv)jv < ||uiv||o,Jv||/||o,JV- Using the 
ellipticity assumption, we obtain the estimates 

||ujv||o,jv < A _1 ||/||o,v 

\\u N \\l N < [(Ae )- 1 +A- 2 ]||/||^. 

By the finite elements theory [TJ, there exists a constant 7 independent of N 
such that for every / € Hi(U d ), ||/||o,v < 7||/||i- Therefore, it is enough to take 
c = 7 max{A- 1 ,(Ae )- 1 + A- 2 }. □ 

Since the operator Ln is the generator of a random walk in T^, the solutions of 
(|B.1|) satisfy the maximum principle: 

inf A -1 f(x) < inf un(x) < inf un(x) < sup A _1 /(a;). 

In particular, for / continuous, ||ujv||oo 
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